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THEORY  OF  SOLID  STATE  CYCLOTRON  MASER 

I . INTRODUCTION 

1 2 

The  cyclotron  maser  radiation  ’ in  an  electron  beam  in  vacuum  originates 
from  the  relativistic  variation  of  electron  mass  (hence  the  cyclotron  fre- 
quency) with  velocity.  In  semiconductors  such  as  InSb,  the  effective  mass 
of  the  electrons  depend  on  velocity  due  to  the  nonparabolic  nature  of  the 

energy  bands.  This  property  can  be  utilized  to  make  a solid  state  cyclotron 
3 

maser  (CM)  in  such  semiconductors.  The  maser  will  operate  in  the  submillimeter 
range  and  provide  a low-power  source  of  radiation  in  this  frequency  range. 

4 

As  shown  by  Kane,  the  energy  (W^)  of  electrons  in  the  conduction  band 
of  InSb  is  related  to  momentum  (p)  by 

2 2vV/2 

"»'L{Es  +"/'!  -\l/2 

* O 

where  mQ  is  the  effective  mass  of  the  electrons  at  the  bottom  of  the  conduction 

C ^ _ 

band  and  E the  band  gap.  From  this  equation  and  the  relation  v =V»W  it 
o p b’ 

follows  that  v ■ p/m^  (1+p  /m^  Vg  ) = p/m  , where  the  parameter  m is 
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• In  InSb  E “0.24  ev.  v = 1.30  X 10  m/sec  and 
g g 

* v 

m Q “ .014  where  me  Is  the  free  electron  mass.  Because  of  small  values  of 
* 

mQ  the  solid  state  cyclotron  maser  can  operate  in  the  submillimeter  frequency 

range  with  relatively  weak  magnetic  fields  of  the  order  of  several  kilogauss. 

If  a thin  InSb  sample  at  77°K  is  placed  in  a magnetic  field  B and 

o 

electrons  are  injected  in  the  conduction  band  at  an  angle  to  the  magnetic 

field,  the  electrons  will  move  in  helical  trajectories.  Initially,  the 

phases  of  the  electrons  in  the  cyclotron  orbits  are  random  and  no  radiation 

is  emitted.  But  phase  bunching  can  occur  due  to  the  dependence  of  cyclotron 
Note:  Manuacript  (ubmitted  November  16,  1978. 
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frequency  (1  - e Bo/m  on  the  electron  energy.  The  electrons  that  are 
decelerated  In  the  wave  electric  field  rotate  faster  and  accumulate  phase  lead 
while  the  electrons  that  are  accelerated  rotate  slower  and  lag  in  phase . This 
results  in  phase  bunching  and  the  electrons  radiate  coherently  at  frequency 
u)  *»  sfl  where  s is  an  integer.  This  mechanism  of  maser  radiation  is  different 
from  the  optically  pumped  "cyclotron  maser"  proposed  earlier.^’ ^ 

Theoretical  calculations  for  cyclotron  maser  in  vacuum  have  been  given 

g 

for  two  different  configurations,  (1)  the  waveguide  structure  and  (2)  the 

9 

cavity  structure.  In  the  first  configuration,  the  electromagnetic  wave  grows 
as  the  result  of  an  instability  driven  by  the  electron  beam.  It  corresponds  to 
travelling  wave  amplification  in  waveguide  structures.  In  the  second  model,  the 
electron  beam  Interacts  with  the  constant  amplitude  standing  wave  of  a cavity 
structure.  It  corresponds  to  beam  sustained  oscillations  in  a finite  Q cavity. 

Recently  a theoretical  treatment  of  cyclotron  maser  in  solids  was  given  by 
3 

Kalmykov  et  al.  for  the  waveguide  structure.  In  their  paper,  the  dispersion 
relation  was  derived  from  the  Boltzmann's  equation  in  which  the  collislonal 
integral  was  disregarded.  Conditions  for  maximum  wave  growth  were  obtained  from 
the  dispersion  relation  and  the  feasibility  of  a solid  state  cyclotron  maser  was 
demonstrated . 

A major  difference  between  the  vacuum  cyclotron  maser  and  the  solid  state 

cyclotron  maser  is  the  effect  of  collisions  in  the  latter.  Collision  in  solids 

will  be  a serious  obstacle  for  the  phase  bunching  needed  for  coherent  readlation. 

The  electrons  will  remain  in  orbit  for  a distance  close  to  the  mean  free  path 

length  X . For  example,  in  InSb  near  77°K,  X is  of  the  order  of  10  - 100  urn.^’*® 
z z 

For  electron  velocity  of  the  order  of  10^m/sec,  the  Larmour  radius  r.  » Vj_/Q  and 

L 

the  pitch  of  the  spiral  X£  = 2tt  v^/Q  are  of  the  order  of  .1  |im  and  1 um,  respective- 
ly,  for  a magnetic  field  — 5 kOe.  vj^  and  v^  are,  respectively,  the  components 


of  electron  velocity  perpendicular  and  parallel  to  the  magnetic  field.  r.  and 

I* 

\c  are  both  much  smaller  than  the  mean  free  path  and  a large  number  of  turns 

(X  /X  ) of  the  electron  spiral  occurs  within  X . Thus,  electron  cyclotron  maser 
z c z 

interaction  can  take  place  If  the  Interaction  length  is  of  the  order  of  X^,  but 
the  conditions  for  cyclotron  maser  interaction  rapidly  deteriorate  when  the 
interaction  length  goes  much  beyond  X^.  In  this  regard,  a solid  state  cyclotron 
maser  in  the  cavity  configuration  is  expected  to  offer  a very  significant  ad- 
vantage over  the  waveguide  configuration  because  the  cavity  configuration  requires 
a much  shorter  interaction  length.  The  reason  is  as  follows.  In  the  waveguide 
case,  the  beam  interacts  with  an  electromagnetic  wave  which  grows  from  the 
noise  or  near  noise  level,  while  in  the  cavity  case  the  beam  interacts  with  a 
large  amplitude  standing  wave  which  has  been  built  up  and  stored  in  the  cavity. 

As  a result,  the  interaction  is  much  stronger  in  the  latter  case,  or  in  other 
words,  the  required  interaction  length  is  much  shorter  in  the  cavity  than  in 
the  waveguide  in  order  for  the  beam  to  lose  the  same  amount  of  energy.  This 
is  reflected  from  the  fact  that  in  vacuum  cyclotron  maser  experiments,  an 
oscillator  (cavity1*  is  generally  shorter  than  an  amplifier  (waveguide)  by  one 
order  of  magnitude. 

Motivated  by  the  above  consideration,  here  we  formulate  a detailed  solid 
state  cyclotron  maser  theory  in  cavity.  The  formalism  and  the  physical  ex- 
pressions to  be  derived  are  considerably  different  from  those  of  the  waveguide 
3 

structure  although  the  basic  mechanism  is  similar.  The  exact  spatial  field 
variation  has  been  incorporated  in  our  calculation  and  the  electron  Larmour 
radius  has  been  kept  arbitrary.  This  allows  us  to  examine  interactions  at  the 
non fundamental  as  well  as  the  fundamental  cyclotron  frequencies.  More  lm- 


3 


portantly,  the  effect  of  collisions  has  also  been  Included  in  our  model.  It 
will  be  shown  that  the  collisional  effect  Is  such  a dominant  limiting  factor 
in  the  solid  state  cyclotron  maser  that  the  short  interaction  length  afforded 
by  the  cavity  structure  could  be  a decisive  advantage  over  the  waveguide  structure. 

In  section  II  we  calculate  the  linear  response  of  an  annular  electron 
beam  as  it  Interacts  with  the  cavity  modes.  The  dynamics  of  the  electron 
beam  is  determined  from  the  Vlasov  equation  while  the  collisions  are  treated 
with  an  approximate  model.  In  section  III  we  calculate  the  electron  beam- 
cavity  mode  coupling  coefficient,  beam  energy  gain  and  the  threshold  beam 
power  necessary  to  sustain  oscillations  on  the  basis  of  a cold  beam  assumption. 

We  also  show  the  dependence  of  these  quantities  on  the  various  parameters  such 
as  beam  position,  beam  energy,  cavity  length  and  the  magnetic  field.  In 
section  IV  we  discuss  the  results  from  our  theory  and  suggest  some  experimental 
configurations  to  observe  the  oscillations. 

J 
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II.  MODEL  AND  FORMULATION 


Fig.l  shows  the  configuration  of  the  electron  cyclotron  maser  system 
under  consideration.  An  annular  beam  of  electrons  is  guided  by  the  magnetic 
field  Bq  along  helical  trajectories  Inside  a circular-crosssection  InSb 
cavity  (radius  R and  length  L < X^).  The  electrons  are  injected  at  an  angle 
to  the  magnetic  field  such  that  a major  part  of  their  kinetic  energy  is  in 
the  form  of  transverse  gyromotion  and  the  rest  in  the  form  of  axial  motion. 

The  axis  of  the  trajectories  is  along  the  cavity  (z-axis).  The  cyclotron 
orbits  may  or  may  not  encircle  the  axis  of  the  cavity,  depending  on  how  the 
beam  is  formed.  In  Fig.l  the  second  type  of  orbits  are  shown.  We  make  the 
following  simplifying  assumptions  to  obtain  the  linear  response  of  the 
electron  beam:  (i)  the  beam  is  sufficiently  weak  so  that  its  self  electro- 
static and  magnetic  fields  are  small  compared  with  the  cavity  fields,  (ii)  the 
cavity  fields  are  of  first  order  with  respect  to  the  applied  magnetic  field 
Bq  and  the  perturbed  electron  distribution  f^  caused  by  these  fields  is 
of  first  order  with  respect  to  the  initial  distribution  function  fQ, 

(iii)  the  distribution  function  and  the  cavity  fields  are  independent  of  the 
azimuthal  angle  6,  and  (iv)  the  electron  collision  frequency  is  much  smaller 
than  the  electron  cyclotron  frequency. 

The  cavity  modes  may  be  classified  as  TE  or  TM  modes.  Cyclotron  maser 

interaction  is  much  stronger  for  TE  modes  than  the  TM  modes. * Thus  we 

consider  only  TE  modes  and  in  accordance  with  assumptions  (1)  and  (ill) 

restrict  our  attention  to  TE  modes  given  by 

onm  ° J 


E^1}  - E-  J.  (k  r)  sin  k z cos  ait  , 
u oin  z 

o 

B^  - (k  /«»)  E„  J (k  r)  cos  k z sin  tut  , 
r z bin  z 

o 

Bz  “ "(kn/u>)  E0  JQ(knr)  sin  kzz  cos  u>t  , 


where 


Hie  electron  distribution  function  f is  determined  by  the  equation 


+ <4> 

m V ® / V / coll 

where  I'd f/3t\  is  a collision  term  whose  form  will  be  specified  later. 

\ hl 

On  the  basis  of  the  two  different  time  scales  assumed  [[assumption  (iv)  ^ , 
we  may  separate  the  distribution  function  f into  a slowly  varying  component  f 
representing  the  collisional  relaxation  of  the  zero  order  distribution  function 
and  a fast  varying  component  representing  the  perturbation  caused  by  the 

wave  field.  A characteristic  of  the  cyclotron  maser  system  is  that  the  cavity 
structure,  rather  than  the  electron  medium,  determines  the  properties  of  the 
wave  (spatial  profile  and  dispersion  relation,  etc.).  Thus,  only  those 
electrons  which  interact  resonantly  with  the  wave  are  of  importance  while  the 
effects  of  nonresonant  electrons  can  be  neglected.  As  will  be  demonstrated  in 
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section  III,  resonant  electrons  occupy  a vary  narrow  region  in  the  velocity 

space.  Hence  the  initial  distribution  function  of  interest  to  the  problem 

would  be  one  in  which  all  the  electrons  fall  in  the  resonant  region. 

Further  it  is  reasonable  to  assume  th-  if  an  electron  suffers  a collision 

in  its  path,  it  will  be  scattered  off  the  resonant  region  and  consequently 

its  dynamics  will  no  longer  be  of  interest.  Thus  we  may  approximately  write 

f as 
s 

fg(r,  P,  t)  = fQ(r,  p)  exp  [ - v(t  - tQ)  J , (5) 

where  t is  the  time  an  electron  first  enters  the  cavity,  v * v /X.  the 

o z ' z 

collision  frequency  and  fQ(r,  p)  any  function  which  satisfies  the  zero 

order  Vlasov  equation^  in  the  absence  of  collisions.  Since  the  cavity  length 

L < X , very  few  collisions  occur  and  the  approximate  form  for  f in  Eq.(5) 
z s 

should  be  adequate  for  our  purpose. 

On  linearizing  Eq.(4)  according  to  the  ordering  scheme  in  assumption  (ii), 
we  obtain  the  following  equation  for  the  perturbed  distribution  function  f^. 


12 

Note  that  in  Eq.(7),  the  BGK  model  of  collisions  has  been  assumed.  Eq.(6) 
may  be  solved  by  the  method  of  characteristics,  namely,  by  Intergrating  it 


7 


along  unperturbed  trajectories  of  the  electrons.  Eq.(6)  thus  reduces  to 


f(1)(r,'p,  t) 


-2* 
u P 


*]• 


^pi  fBCr  '*  p' , t ' ) , 


(8) 


where  the  t'  integration  is  along  the  unperturbed  orbits.  The  primed 

quantities  r ' and  p'  are  treated  as  functions  of  t'  while  r*  and ^ are, 

respectively  the  values  of  r'  and  p 1 at  t'  - t.  The  lower  limit  of  the 

integration  t is  given  by  t - t -z/v  , i,e.  the  time  an  electron  at  axial 
o o z 

position  z and  time  t first  enters  the  cavity.  Substituting  Eq.(5)  in  Eq.(8) 

and  using  the  relation  t - t - z/v  , we  obtain 

o z’ 


^ t 

p,  t)  = exp[-z/\z]  • ^ dt'  e [_^ ( 1 ^ ( r ’ , t') 


t-z/v 


+ xt<1> ]• 


P'  o 


(9) 


In  the  absence  of  collision  f^  is  given  by  the  integral  alone  on  the  right- 
hand  side  of  Eq.  (9) . Coliisions  reduce  f^  by  the  factor  expt-z/X^]. 

The  methods  for  evaluating  the  integral  in  Eq.(9)  are  standard11  and  the 
result  may  be  written  in  the  form 


:(D 


fa>+  f(D 

T “ 


(10) 
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A 


where 


(11) 


with 


0 “ eB/m  - Q/V 
o o 


r . - p,  /m  n , 

L I!  o o * 

X - 1 - exp  ^i(tu  - kzvz  - s 0)  ,/vJ 


G88'(X)  - d J8(X>/dX 


(12) 


In  Eq.(ll),  Imf A ] indicates  the  imaginary  part  of  A.  q>  and  9 are, 
respectively,  the  polar  angles  of  the  momentum  and  position  vectors,  p^  and 
Pz  denote  the  components  of  momentum,  p*  ■ m*  v",  perpendicular  and  parallel, 
respectively  to  the  external  magnetic  field.  ffX)  is  given  by  equation  (11) 
with  te  replaced  by  -u>.  In  obtaining  equation  (11)  we  have  used  the  following 
relation13 


is  0i 

e w 


J-OD 


is'0, 

Js+s,(x2)  Js,(x3)e 


where  x^,  Xj,  x^,  9^  and  9^  ®te  related  through  the  triangle  shown  in  Fig. 2. 
The  perturbed  azimuthal  current  Jg1’)  is  given  by 

J01):*  'e  j f(1)v0d3p  - -e  j P±  dp^  ^ dp2  ^ d£  f(1)  v^  sin  £ (1 

0 -•  0 
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where  CD  - cp  - 9.  From  Eqs . (10) -(13) , we  obtain  after  integration  by  parts 


over  p.  and  p 

1 z 


(1)=  (1)  (1) 

Je  J9f  + V 


(14) 


where 


-e2E 


* 

. 4m  a) 
L o 


9 i(k  z-cut)-zA 
o z z ^ ^ 


» +*  2tt 

( \ (a  , 

. \ p dp  ; dp  j dec  f i 


J do  f i8'  Je-1*8’-1)*  - 
~±  '1  "«  o 

0 -®  0 


} 


| Gss'(knrL)pI  (U)2  - kz  \ )X  . <m  - kzVz)X{2Gss'(kn,rL)+knrLG;s'(knrL)l 


,3  *2  2 2 YTl 

V mo  vg^  ^ 


iz(l-X)(t02  - k2  V2  - k V2T|  /v  '-2 
+ \ z g z g I z 


2*2 
Y m v p 
o g Kz  I 


k 


(15) 


with  Y|  = u)  - kzvz  - sft  and  G^s,(x)  = d Gss,(x)/dx  . is  given  by  Kq.(15) 

with  (u  replaced  by  -oj. 

Equation  (14)  gives  the  perturbed  current  for  a general  distribution 
function  fQ.  In  the  following  we  specialize  to  a particular  distribution 
function  which  is  considered  to  be  the  most  ideal  for  cyclotron  maser  operation.^ 
It  is  constructed  from  the  constants  of  motion  of  the  system  (in  the  absence  of 
collision):  p , p and  L„  « Y m*  v r sin$  - £ B r2.  We  assume  that9 

i Z ff  4 0- 

°) 


£o  ■ k 5 (rL  - 2VeBo 


r2\  °(l-  1°M 

o ) 2np^ 


VPz/ 


(16) 
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which  represents  a cold  monoenergetic  bean  moving  along  helical  trajectories 
with  guiding  centers  distributed  uniformly  on  a cylinder  of  radius  r#  . The 
constant  k is  determined  by  the  normalization  condition 

J fo  2irt  dr  d3p  - N (17) 

where  N is  the  line  density  of  the  electrons  (i.e.  number  of  electrons  per 
unit  length).  Eq.(16)  may  be  rewritten  as 


N «(p„  - P°)  6(P.  - P,°)  I A A a a | 

f - — — 1 -2^—  1/2  |6(cp  - <po)  + 6(cp  - rr  + q>o)  0(r -^XK^-r) 

2nZ  |frZ  - rp(r‘  - O } ' ’ 

(18) 


where 


<p  “ sin 
o 


_1  [f*2  +rr2  ~ 12  r rll  » 


«(x) 


J1 


for  x >_  0 
for  x < 0 , 


*1  “ lr0  ■ V “d  r2  = 'ro+  rL> 


The  Idealized  distribution  function  in  Eq.(18)  leads  to  results  which  can  be 
physically  interpreted  and  provides  valuable  information  for  the  operation  of 
the  solid  state  gyrotron. 


III.  BEAM  POWER  GAIN  AND  THRESHOLD  POWER 

The  time  average  power  gain  for  all  the  electrons  in  the  cavity  is  given 


by 


2n/u>  r L 

P - (o  ^ dt  ] dr  j d«  .T<1)  E<1}  . (19) 

0 0 0 
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i 


• ( 


We  introduce  the  following  notation: 


r 


a 


H (k  r k r )p 

8 n O PL 


s2 


Q (k  r , k r ) 
8 no  n L 


(23) 


(u)t 


kzD  . | m(a')  - r } , 


with 


M(x)  - Cx(l  * e VT  cos  *>  _ VT  e_VT  sin  x]  / (x2  + V2  t1)  f 


N(x)  - Re 


, -vT  ix 
1 - e e 


. -VT  ix 
i e e 

X + ivT  J , 


.(X  + ivT)‘ 

T - 2k  L (1  - e"vT)  / (4k2  L2  + v^2)  . 

z * z 


(24) 


(25) 


(26) 


The  quantities  awi  are,  respectively,  obtained  from  agl  and  a#2  with 
the  following  replacements 

(u  -♦  -a),  A -A'  and  A'  -*  -A  . 

The  functions  Hs(aQ,  aL>  and  Qg(ao,  a L>  are  defined  by 


H * J * (a.  ) I ) > 

S 8 L 8 O L 


' 2W  “l5  + aL  ^s  ‘WV  *L) 

+ 2 aL  Js(aL)  Is-l(ao*  *L)  ' Is+l(ao’  *L)  ^ » 

where 

a2  ^ 2 2 2 1 
Is(ao,aL)  ” ' n ^ da  Jl(a)  a sln*  t (*2  ‘ al)(a2  ‘ *2)  J 
al 

dP 

J§..WC*  V(a)  coa  8’(f-*)  * 


(27) 


(28) 


-1/2 


(29) 
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The  integral  series  in  Eq.(29)  is  evaluated  in  Appendix  A and  found  to 
L 9 


w *l>  • w • 


Eq . (21 ) contains  a simulation  over  all  harmonics  of  the  cyclotron  frequency. 

As  we  shall  see  later,  radiation  is  favored  at  a particular  cyclotron  fre- 
quency near  the  synchronous  condition 

ou-k  v -sO/y  = 0. 
z zo  o o 

Hence  in  Eq.(21)  we  keep  only  one  term.  We  are  mostly  interested  in  the 
fundamental  cyclotron  frequency.  In  later  calculations  parameters  are  chosen 
such  that  s = 1 term  dominates.  In  the  case  v . o,  the  quantities  a assume 


the  form 


■V^l0/Pg^  S^-^s2  /|2  .)  {4  sin2(A/2)  - A sinA} 

+ k L(p2/p2)A^A  sinA  - 2 sin2(A/2)  fj  /a3 

Z R Z I ' 


+ 2Q  . (an  - k L)sin2(A/2)  /A2 


-V  ^X0/Pg^2  [ (U)2  t2-  kz  ^ {2d  + A/A’) Sin2 (A’/ 2) -A  sinA’  } 

+ {A  8inA'  " 2(A/A'>8ln2(A'/2)  | ] / (A2A' ) ^ 

- 2Qs  . (ourr  - k L)8in2(A'/2)  / (AA*)  . 


The  other  two  quantities  and  are  obtained  by  the  replacements 

0)  -*  -u>,  A “*  -A'  and  A1  -*  -A 


in  the  expressions  for  and  o^,  respectively. 


In  Eq.(10)  and  may  be  regarded  as  perturbations  in  the 


electron  distribution  function  due  to  the  forward  and  backward  travelling 


waves  which  make  up  the  standing  cavity  modes.  In  eq.(21)t  the  terms 


:(D 


a . and  a . arise  from  interaction  of  f , with  the  forward  and  the  backward 
si  sr  + 


travelling  waves.  Similarly,  and  terms  are  due  to  the  interaction 


= (1) 


of  f with  the  two  waves. 


In  the  expression  for  a' s,  the  first  term  (proportional  to  Hg  p^g  ) 


arise  from  the  transverse  motion  of  the  electrons  and  the  second  term 


(proportional  to  (ioT  - kzL)Qg)  is  due  to  wave  induced  oscillations.  The 


first  term  is  much  larger  than  the  second  term  unless  p^g  is  too  low.  Hence 


the  coupling  between  the  electron  beam  and  the  cavity  modes  is  essentially 

proportional  to  H (k  r ,k  r ).  The  first  term  may  be  positive  (beam  power 
8 n o n L 

gain)  or  negative  (beam  power  loss)  depending  on  the  phase  factor  A. 


It  is  shown  in  Appendix  5 that  E the  amplitude  of  the  component  of  the 

® 


cavity  electric  field  tangential  to  the  electron  cyclotron  orbit  (see  Flg.l) 
may  be  written  as 


,1/2, 


E,_  =-A  E.  X H (k  r k r ) cos  s ® , 

V ®o  7=0  8 n ° n L 


where  A = 1 for  s = 0 and  A ■ 2 for  s ^ 0.  Thus  E0  H ’ may  be  Interpreted 

o 


as  the  s-th  harmonic  component  of  the  effective  electric  field  in  the  di- 
rection of  the  electron  velocity  and  the  beam  power  gain  is  proportional  to 


V 
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In  the  expression  for  P,  the  guiding  center  position  r q enters  through 

Hs^knro’  knrL^’  Hence  p can  be  “•ximixed  by  maximising  with  respect  to  r 

for  each  combination  of  n and  s.  Numerical  computation  shows  that  the 

optimum  value  of  r b »48R  for  n - s » 1.  This  value  of  r will  be  used  in 
o o 

later  calculations. 

The  total  stored  energy  (W)  in  the  cavity  is  given  by 

W - ttr2L  « J2  (kR)E2  / 4 . (33) 

o 

Let  us  define  a dimensionless  quantity  F 

F = P t/w  , (34) 

as  the  ratio  of  the  total  beam  energy  gain  during  one  transit  time  to  the 
total  stored  field  energy.  From  Equ (21),  (33)  and  (34),  we  get 

F - 0)2  r a/2  Yo  J2  (xn)u>  , (35) 

22  * 1/2 

where  u)p  = (N  e /ttR  « mQ)  and  a - (ctj  + a2  - a3  - 0^) . 

The  quantity  F may  be  positive  or  negative  depending  on  the  value  of  the 
phase  factor  A.  The  beam  generates  electromagnetic  radiation  in  the  range  of 
A for  which  F is  negative. 

An  actual  cavity  has  a finite  Q due  to  dielectric  losses  or  loading  or 
wall  resistivity  or  end- faces  not  being  perfectly  reflecting  and  energy  is 
lost  from  the  cavity  at  the  rate 

Pout  " * W/Q  * <36> 

The  electron  beam  generates  energy  at  the  rate 

Pin"‘FW/T  • 07) 
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Thus  the  threshold  condition  for  cavity  oscillations  is  given  by 


■F  Q i it  T 


The  energy  of  the  electrons  in  the  conduction  band  of  the  crystal  may  be 
written  as 

Wb  “ (Yo  - 1)mo  Pg  c2  • <39) 

and  the  electron  beam  power  inside  the  crystal  is 

P.  - N(V  - l)m*  p2  c2  v . (40) 

b o o g zo 

From  Eqs.  (35),  (38)  and  (40),  we  find  that 

Pb  * Pf  , (41) 

where  the  threshold  beam  power  P^  is  given  by 

Pbh  “ •(u,R/c>2  • { W1)Pg  Pzo  Jo(xn^}  • { 4n  e mo2  c5/2e2  ) • (42) 

The  result  for  cyclotron  maser  in  vacuum  is  recovered  by  setting  6 * e^, 

P =1  and  m « m in  Eq.(42).  Since  P = 4.3X10  3 and  m * 0.014  m in  InSb 
g o e g o e 

at  77°K,  the  threshold  beam  power  for  the  solid  state  cyclotron  maser  is  many 
orders  of  magnitude  lower  than  that  of  the  vacuum  case.  Numerical  examples  for 
F and  P^1  are  given  in  the  next  section. 


IV.  NUMERICAL  RESULTS  AND  DISCUSSION 


The  quantity  F (eq.35)  and  the  product  Q P^  (eq.42)  depend  on  the 

parameters  8,  n,  m,  R,  r , \ , Y , v /v  and  0.  For  numerical  examples  we 

O Z O 1.0  zo 


consider  only  the  fundamental  mode  numbers,  l.e.,  s « m ■ n — 1.  In  this 


case  x * k R = 3.832,  k L = n and  A'  » A + 2n  . As  shown  earlier  r should 
n n z o 

be  -48R  to  maximize  H (x  r /R,  x r./R).  From  equations  (35)  and  (42)  it 
s n o n l 

can  be  seen  that  the  radius  of  the  cavity,  R,  may  be  scaled  out  by  nor- 
malizing length  to  R and  frequency  to  c/R.  Wave  vector,  time  and  velocity, 
are  correspondingly  normalized.  The  natural  dimensionless  quantities 
Yq,  P,  A,  iut , k^L,  Hg  and  Qfl  remain  unchanged.  For  a given  L - L/R,  tv  « 
ctv/R  is  fixed  by  Eq.(3c).  F and  Q are  obtained  numerically  for  various 

values  of  the  parameters  L,  Y , P /P  and  0,  *=  cfl/R  (i.e.  A). 

o zo 

Figs.  3a  and  3b  show  typical  plot  of  F as  a function  of  A. 

F becomes  negative  for  A in  the  range  - . 9tt  < A < 2n  . The  lower  limit 

does  not  change  with  L and  shows  slight  changes  only  with  Y , and  P /B 

o io  zo 

The  maximum  negative  gain,  -F  , increases  with  increase  in  Y , P /P 

® o _[  o ZO 

and  L . The  magnitude  of  F decreases  with  vt  as  expected,  ct^,  O^,  0^ 

and  are  also  plotted  as  function  of  A in  fig.  3.  a ^ differs  only  slightly 

from  a2<  In  the  travelling  wave  structure,  the  gain  is  determined  by  the 

3 8 

term  Ctj.  In  this  case  gain  is  negative  only  ’ for  A > 0.  The  phase 

factor  A corresponding  to  -F  lies  near  + .1  n . A shifts  slightly  towards 
m m m 

lower  values  with  increase  in  L , Y and  p /p 

o j_o  zo 


„th 


In  Fig. 4,  Q is  plotted  as  a function  of  the  electron  energy  in  the  crystal 
*22  — 

Wjj  * (Yq  - l)m0  c Pg  for  cyclotron  frequency  D corresponding  to  A^.  Data 

are  shown  for  P, rt/B  ■ 1.5  and  different  values  of  L and  vt  = L/X  . The 
iu  zo  z 

threshold  power  increases  with  decrease  in  L and  increase  in  vt.  q P*h  is 

-3  ”3  — 

of  the  order  of  10  watts  for  electron  energies  ~ 10  ev.  For  L = .05  and 

vt  « .5,  Q P^h»1.8  mW  at  = 13.2  m ev.  Assumming  Q 91  10C,  P^h  a: 

18  M>W.  As  the-  length  of  the  cavity  is  decreased  to  reduce  the  effect 

of  collisions,  the  lowest  eigenfrequency  Increases.  For  mean  free  path 

X^  — 100  |im  and  L « .05,  the  operating  frequencies  at  vt  = 1.0  and  0.5  are, 

respectively,  367  GHz  and  734  GHz.  The  corresponding  cyclotron  frequencies 


£7  for  maximum  beam  energy  loss  are  357.4  GHz  (Hq  = 1787  Oe)  and  715  GHz  (Hq 
3575  Oe). 


We  make  some  suggestions  for  possible  experiments  to  observe  solid-state 

cyclotron  maser  radiation.  The  cavity  might  be  in  the  form  of  a disc  with 

polished  ends  and  metallized  edges.  The  thickness  of  the  disc  should  be 

less  than  a mean  free  path  of  the  electrons.  The  cavity  is  placed  in  a 

uniform  magnetic  field.  One  end-face  of  the  disc  is  bombarded  by  an  electron 

beam  at  an  angle  to  the  field  such  that  the  electronsget  into  the  conduction 

band  with  appropriate  transverse  and  axial  velocities.  Another  possibility 

is  to  use  the  structure  under  forward  bias  condition  with  a doping  profile 

14  + + 

similar  to  that  of  a Read  diode  (n  pi  p ) such  that  there  is  a very 
narrow  accelerating  region  of  high  d.c.  electric  field  and  a longer  Intrinsic 
region  where  the  d.c.  electric  field  is  negligible.  Cyclotron  maser  radiation 
takes  place  in  this  region. 
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V.  CONCLUSION 

In  this  paper  we  propose  the  operation  of  a solid-state  cyclotron 
maser  in  the  cavity  configuration  and  demonstrate  its  feasibility  with  a theo- 
retical analysis.  We  have  derived  analytical  expressions  for  both  the 
electromagnetic  power  gain  and  the  threshold  beam  power.  The  dependence  of 


these  two  quantities  on  the  various  parameters  such  as  the  length  and 

radius  of  the  cavity,  electron  mean  face  path,  the  magnetic  field,  initial 

electron  beam  energy,  radial  and  axial  eigenmode  number,  cyclotron  mode,  etc. 

has  been  shown.  In  order  to  minimize  the  effects  of  collision,  the  cavity 

length  should  be  less  than  a mean  free  path  of  electrons.  This  sets  a 

lower  limit  to  the  frequency  that  can  be  used.  It  is  seen  from  Eqs. (35)  and 

(40)  that  a better  performance  is  obtained  with  a material  having  large  mean 

free  path  \ , small  band  gap  i.e.  v and  small  effective  mass  m*  . The 
z g o 

cavity  length  can  be  increased  and  the  frequency  of  operation  lowered  if 

is  larger.  Since  the  cyclotron  frequency  is  proportional  to  B /m  , smaller 

o o 

magnetic  fields  are  required  for  materials  with  smaller  mQ  . Furthermore, 

decreases  as  m*2  . A smaller  value  of  v is  also  desirable  because  P*h 
bo  g b 

is  lowered  and  electron  velocity  needed  to  obtain  a given  value  of  V is 

o 

also  reduced.  Our  results  are  based  on  a linear  theory.  A nonlinear  theory 
is  necessary  to  calculate  the  efficiency  and  the  saturation  power  level  of 
the  system.  This  :?ill  be  the  subject  of  a future  investigation. 

The  authors  are  indebted  to  Dr. K.L. Davis  for  many  helpful  discussions. 
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Appendix  A 

EVALUATING  THE  INTEGRAL  SERIES  1,  (o0,  aL) 


The  integral  series  IM  (oq,  aL)  appearing  in  Eqs.  (27)  and  (28)  is  defined  as 


(,  (°0>  a L ) = S da  sin  d>0Jl  — ai  ) (fl2  - «2)  ] *' 


* £ laL)Jt.W  cos  s'  y - *0  . 


a,  -|a0  ~«J. 
«2  - a0  + °L> 


4>0  =sin  1 1 (a7  + al  - ia^)/2aa^] 

Inserting  Eq.  (A.2)  into  Eq.  (A.l)  and  using  the  Bessel  function  identity 


Js  (w)  COS  S*  - 52  COS  Jflr. 


where  »•  — ( u 1 + vl  — 2Mvcosa)'^  and  S'  = cos  *1(m  — vcosc»)/w]  , we  reduce  Eq. 
(A.l)  to  (after  some  algebra) 

s J‘  (o° } f . ^(o)(a2  +a2  -«o)J  I «2  +ao  ~flJl 


11  * * K ° ' Jl  ‘roL  -Oo' 

’ w v,  * Kfl2  -«?)(«!  -a2)]1'2' 


7a0  “L 


To  carry  out  the  integration  in  Eq.  (A.4),  we  replace  the  variable  of  integration  a with  x,  where 
a is  defined  through  the  equation 

a — (a£  + o[  — 2 a0aL  cos  x)  xf2. 

Again  after  some  algebra,  we  obtain 


APPENDIX  A 


^f°l  j dx  J'  1 (°0  + ~ 7a0aL  C05  *)  l/?) 

w 0 ("<?  + aL  ~ 2<3q  "/.  cos  •*) 

• (o^  — a0  cos  x)cos  sx 

1 ' d w 

• ” w ^a0  ) £ rfjrA>  f fro  + "/  ~ 2a0fli  cos  x)  l/2l  cos  «. 

(A3) 

Using  tabulated  integral  formula we  obtain 


Afro-"/)  - A1  ("o) A ("/.)■ 


(A3) 


appendix  b 


L 


The  component  of  the  cavity  electric  field  tangential  to  the  electron 
cyclotron  orbit  is 


E j — Et  cos  a 
” EWJX  (kMr)  cos  o. 


(Bl) 


where  the  z-dependence  of  E+  has  been  neglected.  We  may  express  E6  in  terms  of  r0.  rL,  and 
«5  through  the  following  geometrical  relations  (see  Fig.  1): 


' - ('0  + d ~ 2r0 rL  cos  <t>)  1/2 
COS  a - [rL  + r0  cos  (jr  - <t> ) )/r 
(fL  - r0  cos  0 y ('0  + d - 2 r0rL  cos  <t>)  1/2 


Thus. 


“ EtOJl  ('o  +d  -*rorL  cos  $),/2J(rt  -r0cos«) 
* Id  + r\  - 2rQrL  cos  -,/2 


“ “ £#0fc.  1 37-  2o  l*,  ('o  + d - 2 r0rL  cos  0)  ,/2J. 


(B2) 


Expanding  in  terms  of  the  sinusoidal  harmonics  of  <f>  and  noting  that  E^  is  an  even  func- 
tion of  we  obtain 


E+  - ]£  Et#  cos  s<f> 

9 — 0 


(83) 
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where 

f-j  <**♦««* 

_ W 

QE, o (k„tr  > | *M>  l*«  Oj?  + '2  " 2'o'l  cos 

I 1.  5 -0 

® "1  2.  * ^ 0. 

Using  tabulated  integral  formulae}3  we  obtain  from  Eq.  (B4), 

Eir  “ _ 


- - QEnH^  ttvo-Vi.) 
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Fig.  1 — End  view  of  the  electron  cyclotron  maser.  Guiding  centers  are  uni- 
formly  distributed  on  a circle  of  radius  r0.  'Hie  point  0 is  the  axis  of  sym- 
metry. Hie  circle  of  radius  r^  is  an  arbitrarily  chosen  electron  orbit  and  c 
is  the  guiding  center  of  this  electron.  Eg  the  electric  field  of  the  cavity  and 
E,p  the  component  of  Eg  tangential  to  the  electron  orbit. 


r -i.i 


Fig.  3 (a)  — F (solid  curve)  vs.  A for  n * s = m = 1,  r0*  0.48, 

L = 0.1,  70  - 1.1,  0io//J2o  “ and  VT  = 0.  The  dashed  curves 
marked  c^,  a2»  a3  and  a4  are  plots  of  the  four  components  of 
F as  function  of  A.  Dielectric  constant  e * 17.7  e0. 
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